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Abstract We report on new simulations of the transport of energetic protons 
originating from the decay of energetic neutrons produced in solar flares. Because 
the neutrons are fast-moving but insensitive to the solar wind magnetic field, 
the decay protons are produced over a wide region of space, and they should 
be detectable by current instruments over a broad range of longitudes for many 
hours after a sufficiently large gamma-ray flare. Spacecraft closer to the Sun 
are expected to see orders-of-magnitude higher intensities than those at the 
Earth-Sun distance. The current solar cycle should present an excellent oppor- 
tunity to observe neutron-decay protons with multiple spacecraft over different 
heliographic longitudes and distances from the Sun. 



Keywords: Flares, Energetic Particles; Energetic Particles, Neutrons; Energetic 
Particles, Protons; Energetic Particles, Propagation 



1. Introduction 

Solar flares remain some of the best laboratories in the solar system for study- 
ing particle acceleration, production and interaction. The particles accelerated 
in flares either interact with the solar chromosphere, producing observable X- 
ray and 7-ray emission, or escape into the interplanetary medium. In 1951, 
Biermann, Haxel and Schluter first suggested that solar flares could produce 
energetic neutrons, and evidence of these neutrons was subsequently observed in 
the 2.2-MeV neutron-capture 7-ray line, first by |Chupp et aL|(|1973p and many 
times since. The first direct observations of solar-flare neutrons were made by 



Chupp et al. (1982). Detailed models of neutron production, transport, inter- 



action and decay in flares have been developed (e.g. Lingenfelter et al. [19G5, 



Hua and L ingcnfcltcr , 1987, jMurphy et aT. I2007[) . allowing the number of inter 



acting protons and escaping neutrons from a particular flare to be deduced from 
the 2.2 MeV line fluence. The neutrons which escape the Sun undergo beta decay 
with a mean lifetime at rest of 886 seconds (about 14 minutes): 

n° — > p+ + e" +V e (1) 
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creating neutron-decay protons and electrons in the heliosphere. 

The escaping neutrons and their decay products have been observed by both 
ground-based detectors and space-based instruments near the Earth and around 
the heliosphere ( Chupp et al. 119821 IDebrunner et al\ 11983ft . Since the travel 



time between the Sun and Earth is >20 minutes for particles with energies less 
than 100 MeV, most of the neutrons decay at some time during transit to 1 AU. 
The neutron-decay protons (Evenson, Meyer, and Pyle, 1983||Evenson et al., 1990[)[ 



and also electrons flDroge, Ruffolo, and Klecker, 1996 ) typically arrive before the 



original flare-accelerated protons (protons of the same energy accelerated in 
the flare site rather than originating from neutron decay), and before particles 
accelerated by the shock driven by the associated coronal mass ejection (CME), 
because they make part of their journey as neutrons that travel in a straight line 
and also do not lose energy 

Following the early work by Roelof (1966), several studies have investigated 
the effects that transport in the interplanetary magnetic field (IMF) can have 
on observations of neutron-decay protons, including their time- intensity profiles, 
energy spectra, and pitch-angle distributions (PADs). The first measurements of 
neutron-decay protons found the PADs were nearly isotropic, and, by fitting var- 
ious aspects of the data, obtained scattering mean free-paths of A = 0.3-0.5 AU 



(Evenson, Meyer, and Pyle, 19831, A = 0.2 ± 0.005 AU (Evenson et al, 1983), 



and A ~ 0.3 AU ( Evenson, Kroeger, and Meyer, 1985[ ). In a more complete treat- 



ment Ruffolo (1991) modeled the transport of neutron-decay protons with 25 to 
150 MeV and fit the data from events summarized by |Evenson, Kroeger, and Meyerj 
(|1985ft . He considered pitch- angle scattering and adiabatic focusing of the pro- 
tons and, by fitting "distance-traveled" distributions deduced from data obtained 
by the Third International Sun-Earth Explorer (ISEE-3), obtained the best- fit 
scattering mean free-paths for two events of A = 0.26 ± 0.05 AU and A = 0.37 ± 
0.06 AU. Recently, |Agueda et al. (|2011[) modeled the transport of neutrons and 



their decay products to compare with observations of Solar Cycle 23 events. 

Previous in situ studies of neutrons and their decay products were limited to 
a single vantage point, but new multi-point instrumentation is providing a dif- 
ferent perspective on the geometry of energetic particles in the solar system. The 



canonical view (Reames, Kallenrode, and Stone, 1991 Nitta et al., 20061 holds 



that CME shocks distribute particles widely in longitude, and flare-associated 
protons and electrons cover a narrow (~40°) range in longitude. However, re- 
cent studies have found that small flares can produce particles over ~80° lati- 



tude (Wibberenz and Cane, 2006 Wiedenbeck et al., 20091, while new imaging 



technology is routinely showing previously-undetectablc CMEs associated with 
small flares. Understanding the longitude-distribution of neutron-decay protons 
may help unravel the diverse sources of energetic particles in the heliosphere. 
The twin spacecraft of the Solar TErrestrial RElations Observatory (STEREO, 
IKaiser et aU I2TJTJ8) were launched in 2006 into orbits that lead (STEREO- 
Ahead) and trail (STEREO-Bchind) the Earth. The Advanced Composition 
Explorer ( ACE. [Stone et all 11998ft . Solar and Heliospheric Observatory (SOHO, 
Domingo et al.\ I1983|) . and Wind spacecraft ( Acuna et al., 1995] ) are in orbit 



around Earth's LI Lagrangian point, making multi-point studies covering more 
than 180° in longitude possible in the current solar cycle. 
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2. Description of the Model 

To better understand what these spacecraft should observe in neutron-decay pro- 
tons, we have created a model of neutron and proton transport in the heliosphere 
which determines the trajectories of a large number of individual particles by 
integrating their equations of motion. Neutrons are initially injected as a delta 
function in space and time at the solar surface, for simplicity set at 0° latitude 
and longitude. The energy spectrum and angular distribution of the injected 
neutrons follow two of the examples that were investigated by |Murphy et al.\ 
(|2007p . and illustrated in their Figure 15, where the flare-accelerated proton 
spectrum was assumed to be a power law with a spectral index of -4, and the ac- 
celerated composition was characterizedby He/H = 0.5 along with corresponding 
enrichments in the other species with atomic number Z>2. 

|Murphy et al.\ (|2007[) varied two additional parameters: the mean free-path 
for pitch-angle scattering (PAS) of accelerated protons in the solar atmosphere 
(which we label as A s to distinguish it from the scattering mean free-path in inter- 
planetary space, which wc label as A) and the magnetic convergence parameter 
(<5). In |Murphy et at. (|2007|) A s is a dimensionless ratio of a) the distance for 



an arbitrary initial angular distribution to isotropizc to b) the half-length of 
the loop on which the protons are trapped. They considered cases of no PAS 
(A s — > oo), intermediate PAS (A s = 300), and high scattering (A s = 20). The 
magnetic convergence parameter characterizes the magnetic field strength below 
the transition region, which is assumed to be proportional to the pressure to the 
power 8. For (5 = there is no mirroring, and the angular distribution of the 
interacting protons is downward-isotropic. For 8 = 0.20 the protons will mirror 
as the field increases with depth, leading to an angular distribution of escaping 
neutrons that is somewhat peaked at 90° to the zenith, especially at energies 
>10 MeV. 

In this paper we consider two examples: Case 1 (where A s = 300), and Case 
2 (with \ s — y oo), both with 8 = 0.2. The corresponding spectra of escaping 
neutrons are shown in Figure 15 of |Murphy et al. (|2007p in three angular ranges. 



Wc have approximated the escaping neutron spectra with piece- wise power laws 
covering three energy intervals and three angular ranges from cos 6* = 0° to 
cos 8 = 90° as specified in Tables 1 and 2. Here 8 is the angle relative to the 
surface normal, with 8 = 0° being the direction away from the Sun. We have 
considered neutrons with initial energies between 0.1 and 100 MeV, with flat 
spectra between 0.1 and 1 MeV based on a projection of the curves in Figure 15 of 
|Murphy et~al\ (|2"0"0"T|) . Because our code does not include relativistic corrections 
to the energy loss during transport, wc do not consider particles with energies 
greater than 100 MeV. 

Calculations based on observed 2.2-MeV line fluences have shown that in large 
flares between ~3xl0 28 and ^10 33 protons with energies >30 MeV are acceler- 



ated (Shih, Lin, and Smith, 2009), and the calculations of neutron production 
described above show that for the Case 1 spectra about 0.0054 escaping neutrons 
are produced for each accelerated >30 MeV proton ( jMurphy et al, 20"07 ). For 



the Case 2 spectra where X s — s> oo about 0.0073 escaping neutrons arc produced 
for each accelerated proton >30 MeV. The yields in Tables 1 and 2 are normalized 
to 1 flare-accelerated proton >30 MeV. 
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Table 1. Representation of the Escaping Neutron Spectra: Case l 1 



Energy 


cos 8 = 


0.0-0.3 


cos 8 = 


0.3-0.7 


cos (9 = 


0.7-1.0 


cos e > o 


(MeV) 


Intercept 2 


Slope 2 


Intercept 


Slope 


Intercept 


Slope 


N/P30 3 


0.1-1 


4.51E-5 


0.00 


6.85E-5 


0.00 


7.05E-5 


0.00 


3.51E-4 


1-10 


4.67E-5 


-0.223 


7.83E-5 


-0.429 


9.77E-5 


-0.830 


2.02E-3 


10-100 


5.10E-4 


-1.16 


6.13E-4 


-1.30 


8.95E-4 


-1.68 


2.99E-3 


0.1-100 
















N/P30 3 


1.94E-3 


2.38E-3 


1.04E-3 


5.36E-3 



Case 1 is based on Figure 15b of (Mu rphy et al, 2007| > with X s = 300 and 6 = 0.20 (see text) 
2 The Intercept and Slope entries define power-law representations of neutron energy spectra 
for the indicated energy and angular ranges. 

3 N/P30 is the number of neutrons leaving the Sun for each accelerated proton >30 MeV. 
N/P30 values in the last column and row are integrated over angle and energy respectively. 



Table 2. Representation of the Escaping Neutron Spectra: Case 2 1 



Energy 


cos 8 = 


0.0-0.3 


cos 8 = 


0.3-0.7 


cos 8 = 


0.7-1.0 


cos 8 > 


(MeV) 


Intercept 2 


Slope 2 


Intercept 


Slope 


Intercept 


Slope 


N/P30 3 


0.1-1 


5.42E-5 


0.00 


6.78E-5 


0.00 


8.33E-5 


0.00 


3.86E-4 


1-10 


5.85E-5 


-0.112 


8.46E-5 


-0.422 


9.00E-5 


-0.843 


2.33E-3 


10-100 


7.34E-4 


-1.10 


6.58E-4 


-1.22 


7.34E-4 


-1.52 


4.58E-3 


0.1-100 
















N/P30 3 


3.17E-3 


2.95E-3 


1.17E-3 


7.29E-3 



-"^Case 2 is based on Figure 15a of (Murphy et al, 2007) with A s — ¥ oo and <5 = 0.20 (see text) 
2 The Intercept and Slope entries define power- law representations of neutron energy spectra 
for the indicated energy and angular ranges. 

3 N/P30 is the number of neutrons leaving the Sun for each accelerated proton >30 MeV. 
N/P30 values in the last column and row are integrated over angle and energy respectively. 



At each timestep in the simulation, the escaping neutrons are tested for decay 
by comparing a random number between and 1 with a decay probability given 
by Pa = At/rd where Td is the mean lifetime (corrected for time dilation). Once 
a neutron decays, the trajectory of the resulting proton is followed as it begins to 
feel the effects of the interplanetary magnetic field, assumed to be a Parker spiral 
for a solar wind speed of 400 km s -1 . Rather than computing the gyromotion 
directly, we treat the particles adiabatically, that is, we assume the perpendic- 
ular kinetic energy over the magnetic field magnitude (W±/B) is constant in 
the frame corotating with the Sun. This scheme produces results in excellent 
agreement with schemes where the gyromotion is integrated directly, but results 
in a significant savings of computational resources. The proton trajectory is 
followed until it is lost to the Sun, passes the outer boundary at 10 astronomical 
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units (AU, 1 AU is the Earth-Sun distance), or until it has been in transit for 
10 days. 

Our code includes the effects of pitch-angle scattering and adiabatic cooling of 
the protons in the expanding solar wind. The particle scattering is included in a 
manner similar to the particle decay, where a scatter probability is compared to 
a random number for each timestep. The form of the scattering mean free-path is 
derived from the quasilinear theory ( |Jokipii, 1966| |Giacalone and Jokipii, 19 99 ) : 

W,.)-ai£) J (^) J , 

where r is radial distance, v is particle velocity, B is the strength of the IMF, 
and v is the velocity of a 1-McV proton. 

Our value for the scattering mean free-path A ranges from 0.17 to 0.22 AU for 
25-150 MeV protons at 1 AU and gradually decreases closer to the Sun, values 
that agree well with the previously observed PADs in neutron-decay proton 
events. The recent study of neutron-decay protons by Agucda et al. (2011) 
assumed scatter-free propagation, in which the decay-protons quickly focus into 
a beam aligned with the IMF. In such a beam, protons cross a detector at 1 
AU only once. In our simulation, 10-MeV protons cross 1 AU an average of 6-7 
times. 

In an earlier paper reporting on simulations of solar energetic particle (SEP) 
transport, |Chollet, Giacalone, and Mewaldt|(|2010p found that the proton energy 
deposition at 1 AU was constant to within a factor of two as the scattering mean 
free-path was varied from A = 0.01 to A = 1 AU. It was also relatively insensitive 
to variations in the radial dependence of A. This is because when pitch-angle 
scattering is increased the protons cross 1 AU more times on average, but they 



also lose more energy (see Figure 7 in Chollet, Giacalone, and Mewaldt (20 1U)) ). 



If A is increased, we expect about the same proton fluence, but distributed over 
a shorter time period, with higher average intensity. As a result, we believe 
the results of our simulations are relatively insensitive to the details of our 
description of pitch-angle scattering. The events observed to date suggest that 
the inclusion of some scattering scheme is necessary if the simulation is to match 
the details of real events. 

Our simulation outputs the times particles cross a sphere at a given distance 
from the Sun along with their positions and energies, allowing intensities and 
fluences to be calculated at the simulated spacecraft. Calculations based on ob- 
served 2.2-MeV line fluences have shown that in large flares between ^3xl0 28 and 
^10 33 protons with energies >30 MeV are accelerated ( |Shih, Lin, and Smith, 2009| ).| 

The calculations of neutron production described above show that for the 
Case 1 spectra about 0.0054 escaping neutrons are produced for each accelerated 
>30 MeV proton ( |Murphy et al, 2007 1 . For the Case 2 spectra where \ s — > oo 



about 0.0073 escaping neutrons are produced for accelerated proton >30 MeV. 
Therefore, we have normalized the intensities and fluences to the Case 1 and 
Case 2 yields of 5.36 x 10 29 and 7.29 x 10 29 injected neutrons expected for 10 32 
accelerated protons >30 MeV, somewhat greater than the median value observed 
in real flares. [At least 23% of the Ramaty High Energy Solar Spectroscopic 
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Figure 1. (a) Intensities of 1 to 10 MeV protons at 1 AU for spacecraft at different longitudes 
following a central-meridian 0°-latitude 7-ray flare with 10 32 accelerated protons >30 MeV 
and Case-1 characteristics (X a = 300 and 8 = 0.20; see Murphy et al. (2007)). The time T 
= marks the time the flare photons would be observed at Earth. To the left of the black 
dotted line, the intensities are binned in 10-minute increments, and to the right of the dotted 
line the intensities are binned in 2-hour increments. For all but the orange curve the mean 
statistical uncertainty is ~18% for the 10-minute points and ~8% for the 2-hour points. For 
the orange curve the mean statistical error is ~43% for the 10-minute points and ~19% for 
the 2-hour points. The different colors and symbols represent the different virtual spacecraft, 
with a legend at the top left indicating the relative longitudes of the spacecraft. The dotted 
line in the legend is the radial direction from the injection point. The sector represented by 
the red x has intensities a factor of 10 lower than the adjacent sector, and so is off the scale in 
this plot, (b) Same as panel (a) but for Case 2 neutron spectra characterized by \ s — > 00 and 
S = 0.20. The mean statistical uncertainties are similar to those in panel (a). Both simulations 
used the same sequence of random numbers. 



Imager (RHESSI) 7-ray flares in which the 2.2 MeV line was detected were this 
large flShih, Lin, and Smith, 2009 1; in some events RHESSI missed some of the 
emission because of shadowing by Earth.] Our total Case 2 yield of 1-100 MeV 
escaping neutrons is about a factor of ^3.5 greater than for similar neutron 
spectra in IHua et ai\ (|2002[) . which was the basis of the transport study by 
|Agueda et aZ.| (|201ip . This is likely because Hua et al. assumed an accelerated 
composition less rich in He and heavier ions, as well as a harder proton spectrum 
(spectral index of -3.5). 



3. Results and Discussion 

Figure 1 shows Case 1 and Case 2 intensities of energetic protons vs. time for a 
simulated spacecraft at 1 AU located at different heliographic longitudes. Each 
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of the seven simulated spacecraft integrate over 45° in longitude and 15° in 
latitude. The 1-AU intensities reflect the greater yield of escaping neutrons in 
Case 2, in most directions the intensities are somewhat higher in Figure lb 
than in Figure la. Increasing the amount of PAS decreases the yield of escaping 
neutrons because more accelerated protons are scattered into the loss cone to 
the chromosphere, so that they interact deeper in the atmosphere while typically 
headed downward. Since the secondary neutrons generally go forward at high 
energy, fewer neutrons escape. When PAS is ignored, more neutrons interact 
while mirroring and so there is increased neutron escape near the solar horizon 
(0 = 90°). 

The intensities in Figure 1 are above the background levels of some currently- 
available instruments over more than half of the longitudinal circle for most 
of a day. The Low Energy Telescope (LET) and High Energy Telescope (HET) 
sensors on STEREO, which observe protons between 2 and 100 McV, have a 
galactic cosmic ray background level of a few times 10~ 5 (cm 2 -s-sr-MeV) _1 
at solar maximum (see also IGomez-Herrero et ail 120101 ILeske et ail 120081 and 
IMewaldt et al.\ 120091) . With this instrument background level, the simulated 
event would be detectable at five of the virtual spacecraft in this simulation 
(225°) for a few hours to a day. The estimated decay-proton intensities for the 



28 October 2003 "Halloween" event would be ~10 times greater. Agueda et al. 
([20TT|) conducted a survey of 1 to 11 MeV Wind-3DP data for eight X-class 
flares, and found no evidence for detection of neutron-decay protons. However, 
their background levels were ^0.04-0.1 cm" 2 sr~ - 1 MeV -1 for 3-6 McV protons 
and ~0.02 cm~ 2 sr _1 MeV -1 for 6-11 MeV protons, and they note that their 
simulated neutron-decay proton spectrum is well below the observed background 
proton spectrum for the only one of these events in which RHESSI observed 
2.2-MeV 7-ray emission. 

Though the intensities predicted by our simulation are above instrument de- 
tection thresholds, for magnetically well-connected events it is difficult to distin- 
guish neutron-decay protons from shock and flare accelerated protons ( jRoelof, 1966[ ).[ 
Events originating in the eastern hemisphere are more likely to be detected 
because the directly-accelerated protons often arrive several hours after the 
neutron-decay products. However, many of the largest SEP events are preceded 
by earlier events that elevate interplanetary SEP levels, which may explain why 
so few detections of neutron-decay events have been reported. 

In Figure 2 we show snapshots of particle positions in the ecliptic plane 
20 minutes and one hour after the flare, with the angles and symbols of the 
virtual spacecraft superimposed. This simulation shows that solar flares provide 
a means of quickly distributing neutron-decay protons (and electrons) over more 
than 180° in longitude within the first hour following the flare. The solar sur- 
face creates a solid horizon that neutrons cannot cross. If the neutron-decay 
protons travel without scattering along interplanetary magnetic field lines, this 
horizon will create a shadowed region in interplanetary space as discussed by 



Evenson, Meyer, and Pyle| (I1983j) and |Agueda ef~aL\ (|2011|) . The scattering in- 



cluded in this simulation allows protons to move against the solar wind flow 
into this shadowed region, albeit with difficulty As a result, a spacecraft like 
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20 minutes 




Figure 2. Snapshots of particle positions at 20 and 60 minutes after the neutrons are injected 
(in 20 minutes a 100 MeV neutron travels just over 1 AU). The left half shows the positions 
of the neutrons, while the right half shows the positions of neutron-decay protons. A legend 
representing the sectors of the virtual spacecraft is superimposed. The neutrons travel in 
straight lines, while the protons spiral while following the magnetic field. This simulation used 
the Case-1 neutron spectra that includes PAS of the accelerated protons with A s = 300 and S 
= 0.20 (see text and Murphy et at (2007)). 

Messenger that happens to be inside the shadow might still be able to detect 
neutron-decay protons that diffuse into the shadow region. 

|Murphy et d] (|2"0UT|) found that more neutrons are produced close to the 
neutron horizon than in the radial direction (see Tables 1 and 2, where the 
minimum neutron intensities are roughly normal to the flare site), and the inter- 
planetary results reflect this. Under the influence of the interplanetary magnetic 
field neutron-decay protons move clockwise from the sector marked with a square 
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Energy (MeV) 

Figure 3. Proton fluence spectra integrated over the first four hours following neutron release 
(assumed to be the peak of the flare) arc shown for each of the virtual spacecraft following a 
central-meridian 7-ray flare with 10 32 accelerated protons >30 MeV. The relative intensities 
at the different spacecraft reflect the time-intensity profiles in Figure 1. Low-energy protons 
below 0.5 MeV have not yet reached f AU. Panel (a) is based on Case f neutron spectra with 
A s = 300 and <5 = 0.20. Panel (b) is based on Case 2 neutron spectra with X B — ¥ 00 and 8 = 
0.20. 

(W67.5°-W112.5°) to the sector marked with an asterisk (W22.5°-W67.5°), and 
from one half of the sector marked with a plus (W45°-W67.5°) to the other 
(W22.5°-W45°), making the plus and asterisk sectors the highest intensities 
early in the event (see also Figure 1). Later in the event, the particles have moved 
even more clockwise into the diamond (E112.5°-E157.5°) and circle (E22.5°- 
W22.5°) sectors, while the square (W67.5°-W112.5°) sector quickly becomes 
depleted. 

Proton fluence spectra at 1 AU integrated over the first four hours following 
the neutron injection are presented for each sector in Figure 3. Most of the low- 
speed protons (with energies less than a few MeV) have not yet reached 1 AU at 
this time. The higher-speed protons reflect the injection spectra and transport 
effects. As mentioned above, |Agueda et a/.|(|20lT|) have simulated neutron-decay 
proton intensities at f-AU using somewhat different neutron spectra and in- 
terplanetary transport conditions. Their intensity spectra are similar in shape 
to our fluence spectra (Figure 3), but our spectra show evidence for adiabatic 
energy loss below ~10 MeV. In addition, their predicted peak intensities are 
greater than ours, presumably because, without pitch-angle scattering, all of the 
protons of a given energy arrive within a narrow time interval. 
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Time (Hours) Time (Hours) 



Figure 4. Proton intensities are shown for spacecraft at 0.1, 0.3, 0.5 and 1 AU following a 
central-meridian 7-ray flare with 10 32 accelerated protons >30 MeV. Panels (a) and (c) are 
based on Case 1 neutron spectra while Panel (b) and (d) are based on Case 2 neutron spectra. 
To the left of the black dotted line, the intensities are binned in 10-minute increments, and 
to the right of the dotted line the intensities are binned in 2-hour increments. The virtual 
spacecraft in this case are in the magnetically- well-connected sector spanning 45° , so the black 
solid lines in this plot are the same as the cyan triangle lines in Figure 1. The mean statistical 
error on the 10-minute points increases closer to the Sun, and is ~10% on the 0.5 AU line, 
~20% on the 0.3 AU line and ~40% on the 0.1 AU line. 



As expected, these simulation results show that neutron-decay protons should 
be more easily detectable at heliocentric distances less than 1 AU. The simulated 
spacecraft in Figure 4 are all in the magnetically-connected sector (45° clockwise 
from radial) but at different distances from the Sun, so the 1 AU results in the 
top portion of Figure 4 are the same as the results for that sector in Figure 1. 
Within 0.5 AU heliocentric distance, neutron-decay protons peak within the first 
hour, and remain well above nominal background levels for a full day. The peak 
fluxes at 0.1 AU, 0.3 AU and 0.5 AU are respectively -1000, -100 and -10 
times higher than at 1 AU. The Solar Orbiter and Solar Probe Plus missions, 
which will approach to within 62 and 9.5 solar radii of the Sun (respectively), will 
carry instrumentation that should detect neutron-decay protons and electrons 
in a number of events. 

Close to the Sun, a shock from the associated CME can overtake neutron- 
decay protons with energies of a few McV or less if their radial transport is 
stalled by pitch-angle scattering (see Figure 2). Thus neutron-decay protons 
may provide seed particles for acceleration by CME-driven shocks, as suggested 
bv lFeldman et al.H|2UTP|) . 
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4. Conclusions 

The current solar cycle should present an excellent opportunity to observe neu- 
tron-decay protons with multiple spacecraft over different hcliographic longi- 
tudes and distances from the Sun. Our simulation suggests that these neutron- 
decay protons are very rapidly distributed over more than 180° in longitude and 
should be detectable by current instrumentation at more than one longitude fol- 
lowing large 7-ray events, provided the pre-event energetic particle background is 
sufficiently quiet. Spacecraft closer to the Sun such as MESSENGER, and even- 
tually Solar Probe Plus and Solar Orbiter, will observe even higher intensities 
of neutron-decay protons. While the decay protons and electrons with energies 
greater than a few MeV quickly escape into the heliosphere, lower-energy protons 
linger close to the Sun, where they are available to be accelerated by a CME 
shock into high-energy particles. 
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